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A Remark on the Kramers Problem 

Alex A. Samoletov^ 



We present new point of view on the old problem, the Kramers problem. The 
passage from the Fokker-Planck equation to the Smoluchowski equation, includ- 
ing corrections to the Smoluchowski current, is treated through an asymptotic 
expansion of the solution of the stochastic dynamical equations. The case of an 
extremely weak force of friction is also discussed. 
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1. Evolution of a physical system can be ordered in niulti-time- 
scales. Details of evolution on short-time-scale do not need for description 
in a closed form of a system evolution on long-time-scale and appears on 
this scale only in an average form. The prototype of such kind physical sys- 
tems is dissipative Brownian motion of a particle in an external potential 
field. In this problem, with the exception of extremely short characteristic 
time scales of random forces, there are two time scales: (1) time scales of a 
particle motion in an external field; (2) time scales of relaxation (rate of dis- 
sipation) of Brownian particle in a media. It is intuitively absolutely clear 
that, if the friction force is strong enough (time of free motion is extremely 
short), then probability distribution of a particle velocity to be rapidly 
relaxed to the Maxwell distribution and on this background a particle po- 
sition will be undergoing to slow process of diffusion. In the following we 
deal with the consideration of approximate reduction of the Fokker-Planck 
equation for phase-space probability density to the Smoluchowski equation 
which deals with probability density of a particle position only. In the op- 
posite case of extremely weak force of friction we have energy (or action) 
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variable as evident slow one. 

2. The Kramers problem consist in mathematical description of ap- 
proximate reductions of the Fokkcr-Planck equation for dissipativc Brow- 
nian motion of a particle in an external field to the two limiting cases: (1) 
to the Smoluchowski equation (extremely strong force of friction) or (2) to 
equation for probability density of energy (or action) variable (extremely 
weak force of friction) These reduction procedures are the prototypes 
of all adiabatic elimination procedures or procedures of separation on slow 
and fast variables 

The Kramers model firstly formulated for kinetics of chemical re- 
actions, consists of a particle of mass m moving in an one-dimensional 
potential field U{x) imdcr influence of a random force f{t) and a linear 
friction force with a constant dissipation rate A. The corresponding set of 
Langevin equations has the following form 

X = u, mil = ~U' {x) — \mu + f{t), (1) 

where the random force f{t) is generalized Gaussian 5-correlated stochas- 
tic processes (white noise) with the following properties (including the 
fluctuation-dissipation relation) 

{fit)) = 0, {f{t)f{t')) = 2XmkBTS{t - 1'); (2) 

(• • • ) denotes average over all realizations of random force. 

The Langevin (l)-(2) dynamics is stochastically equivalent to the 
Fokker-Planck equation for the rate of change of probability density P{u, x; t) 
which has the form (e.g.^^^) 

dtP{u, x; t) = -ud^P + -U'{x)duP + -du[uP + ^duP], (3) 
m mm 

where dt = d/dt, dx = d/dx, and 9„ = d/du. 

Keeping in mind (l)-(3) we formulate the problem in the following 
manner: in the case of extremely strong force of friction beginning with 
(3) or equivalently (l)-(2) to derive the approximate reduction to an equa- 
tion for the rate of change of probability density P{x;t) in the form of 
asymptotic expansion by the parameter A~^: 

dtP{x;t) = -dx[X-'Js + o{X-% (4) 

where Jg is the Smoluchowski current, 

Js = -[U'{x)P{x;t) + kBTdxP{x;t)]/m (5) 

In other words it is asymptotics of strong force of friction on time scales 
Xt » 1. 
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This problem has long history starting since 1940, the date of pub- 
lication of the Kramers famous work^^^. For relevant references including 
reviews of the problem see refs. 2-5. The first treatment of the problem 
has been down in rcf. 6 and the first correct solution has been down in 
ref. 7 and then in refs. 8-10. The works'-^""^'*'' are of importance for the 
following in respect of treatment of the corrections older then which 
break the Fokkcr-Planck structure of (4). Most general treatment of the 
problem has been down in ref. 15. 

3. All of the cited works deal with (3) as the input equation for a 
solution of the problem. The purpose of this paper is to take notice of 
the fact that (l)-(2) are indeed convenient input equations for an answer 
to the problem. With respect to solutions of (l)-(2) we use the method of 
asymptotic expansion by the parameter A~^. In the way, the Fokker-Planck 
type equations (4) to be derived from approximate stochastic dynamical 
equations in each order of A""'^. Moreover, Fokker-Planck equation is an ap- 
proximate equation and in any case must be derived from input dynamical 
equations. Convenient and powerful method of derivation, in particular, of 
the Fokkcr-Planck type equations immediately from stochastic equations 
has been initiated by Novikov^^^^ and then it has been sufficiently devel- 
oped by Klyatskin^^^\ In the following we use this method systematically. 
In this connection it should be pointed out that the Klyatskin-Novikov the- 
ory, generally, interprets a stochastic differential equation in the sense of 
Stratonovich. However, in the case under consideration it does not matter. 
We refer the reader to refs. 16 and 17 for further information. 

The probability density P{x; t) can be written in the form*-^^) 

P{x;t) = {5{x-x(t))), 

where x{t) is a stochastic process and 5{- ■ ■) is 5-function. Differentiating 
this definition by time we obtain the equation 

dtP{x;t) = -d^{x{t)S{x - x{t))) = -d^J{x;t), (6) 

which has the form of a conservation law and is the proforma for an ccjuation 
of the type (4). If x{t) is defined by (1), our problem is in calculation of 
asymptotic expansion of x{t) by A~^ and then the corresponding average 
in (6). Further insight is gained by making the following construction. 
Rewrite (1) for Xt^ 1 in the form 

x{t) = -^A-^U'{x) + ^m, (7) 
where operator A has the form K = d/dt + X, and the Ornstein-Uhlenbeck 
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process is introduced: 

at) = A-V(i) = exp(-A<) f eM^t')f{t')dt'. 

Jo 

Formal expansion of by has the form 

A-if/'(x) = \y ^ ( f-) U'{x) + .-.. 
^ ' A ^ A" \dt^ J ^ ' 

n=0 ^ ' 

Hence, (7) can be written in the form 

X{t)r 



m\ m 



[U"{x)x{t)] 



[U"'{x){x{t)f + U"{x)x{t)] + 



The x{t), x{t), and so on, can be exchided from right hand side of last 
equation repeatedly using iterations of this equation and its time deriva- 
tives. Then, and it is important, we are in need of expansion by A~^ of 
the stochastic process (_{t) or, more precisely, of expansion by A"'^ of the 
average in (6) which involves S,{t). First of all we must remark that the 
derivatives f{t) and so on, have sense only as derivatives of the generalized 
stochastic process /(t)'-^*-' and break the simple Fokkcr-Planck structure 
of (6) as of an second order partial differential equation. It is evident in 
the frame of Klyatskin-Novikov theory Namely in the process of cal- 
culation of the corresponding averages according to rcf. 17 we easy detect 
a complex form of (6) including memory as well as an integral-operator 
structure. Further, in respect of X^^f{t) it is easy to verify'^"^) that cor- 
responding averages in (6) have factor A~^ because the noise (2) intensity 
has the order A. 

Taking into account what has been outlined above the first terms 
of expansion of x(t) by A^^ that lead to the current J(x:t) expansion up 
to order A^'^ (the maximum-order of saving of the simple Fokker-Planck 
structure of (6)) can be written in the form 



x{t) 



1 + 



mA2 



U'\x) 



Substituting last expression into the current J{x;t) (6) and performing 
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averaging exactly follow Klyatskin-Novikov theory^^^^ we obtain 

d{x - x{t))) + oiX-"") 



= (l + ^U"{x)^ Js{x; t) + o(A-3). (8) 

where Js is the Smoluchowski current (5). Asymptotic expansions of x{t) 
and J(x; t)(8) together with (2) lead to conventional conclusion: the Smolu- 
chowski equation is valid if: At > 1, l \U'{x)\ < ksT, l'^\U"{x)\ < ksT, 
- where a length scale I = \fk^Tjrn}? is introduced. (8) contains lowest 
order correction to the Smoluchowski equation*^**^^"^. Higher order cor- 
rections in A~^, involving in averaging time derivatives of the generalized 
stochastic process /(t), lead to break of simple structure of the Smolu- 
chowski equation as a second order partial differential equation. It was 
pointed out also in traditional approach^^""^^'^). 

4. Consider now the case of extremely weak force of friction. This 
case more complicated then previous but not so interesting in calculation. 

Let m = 1 in (1). In the case of extremely weak force of friction and on 
the time-scale At <C 1 the energy E = + U{x) of unperturbed system 
is evident candidate for slow variable. But previously E must be averaged 
over period of relatively rapid dynamical oscillations. It is more convenient, 
however , to consider the action variable J instead of E, J = J{E)M^ Let 
J is action variable averaged over period of rapid dynamical oscillations. 
Then an equation for the rate of change of probability density P{J; t) can 
be written in the form (see (6)) 

dtP{J; t) = -dj{j{t)S{J - Jm, P{J; t) = {S{J - J(t))), (9) 

In usual way(^) and taking into account the change of time-scale of the 
white noise^^'^^ we obtain the equations of motion for slow variables 

j{t) = -\J+-.f{t), V=-X-J + fit), (10) 

where iv = to (J) = dE/dJ is frequency and V is velocity averaged over 
period of dynamical motion. Substituting (10) in (9) we obtain 



dtP{J;t) = -dj 



-xjp+-l-{f{t)vmj-jm 



(11) 



For calculation of the average in right hand side of (11) we can use the 
Klyatskin-Novikov procedure again. Using (9) and the causality condition^^^^ 
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we obtain for functional derivatives 



SVit) 



SJ{t) V 



6 fit) ' Sf{t) u- 

Taking also into account that V^/w = J if J=const, according to ref. 17 
we finally obtain 



dtP{J;t) = dj 



XJP - ^(?3&<5(J - At)) - V{t)dj5{J - J{t))^'^^'^ ' 



uj{J) ' 6 fit) 



Sfit)' 



= dj 



XJP-^P+^djiJP) 



;iJ) 



;(J) 



XJ + XksT 



J 



coiJ) 



dj 



PiJ;t). 



It is exactly the Kramers equation. We can derive corrections to this equa- 
tion but it is slightly more difficult task then in the case of extremely strong 
force of friction and does not take special interest in the context of this pa- 
per. 



5. In conclusion, wc have presented in a simplest framework a unique 
approach to the kinetic equations for slow variables by taking stochastic 
dynamical equations as the input instead of the Fokker-Planck equation. 
We hope that this approach is general enough. 
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